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Abstract 



We study locally self-similar processes (LSSPs) in Silverman's sense. By deriving the 
minimum mean-square optimal kernel within Cohen's class counterpart of time-frequency 
i-^h 1 representations, we obtain an optimal estimation for the scale invariant Wigner spectrum 

d . (SIWS) of Gaussian LSSPs. The class of estimators is completely characterized in terms of 

kernels, so the optimal kernel minimizes the mean-square error of the estimation. We obtain 
the SIWS estimation for two cases: global and local, where in the local case, the kernel is 
allowed to vary with time and frequency. We also introduce two generalizations of LSSPs: 
the locally self-similar chrip process and the multicomponent locally self-similar process, and 
obtain their optimal kernels. Finally, the performance and accuracy of the estimation is 
(T) . studied via simulation. 
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F 1 Introduction 

Scale invariance (or self-similarity), once acknowledged as an important feature, has often been 
used as a fundamental property to interpret many natural and man-made phenomena [5]. Self- 
similar processes have been used successfully to model data exhibiting long memory and arising 
in a wide variety of fields, ranging from physics (turbulence, hydrology, solid-state physics, • • • ) 
or biology (DNA sequences, heat rate variability, auditory nerves spike trains, • • • ), to human- 
operated systems (telecommunications network traffic, image processing, pattern recognition, 
finance, • • • ) |17j . Because they may correspond to non-standard situations in signal processing 
or time-series analysis (non-stationarity, long range dependence, • • • ), scale invariant processes 
raise challenging problems in terms of analysis, synthesis, and processing (filtering, prediction, 
• • • ), and a number of specific tools have, however, developed over the years [T7]. Even though 
there is no single definition of scale invariance, it is often described as a symmetry of the system 
relatively to a transformation of a scale, that is mainly a dilation or a contraction (up to some 
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renormalization) of the system parameters [5]. In mathematical expression, a process is scale 
invariant (or self-similar, denoted 'H-ss') [5], if for all A > 0, 

{(D HtX X)(t) := \~ H X(\t),t > 0} = {X(t),t > 0} (1) 

(where = means equal in all finite dimensional distributions). The index H characterizes the 
self-similar behavior of the process, and a very large variety of methods has been proposed in 
the literature for estimating it [2J, [3J, [TTj . 

Sometimes it happens that the H-ss model is not quite adequate for real world phenomenon, 
and it would be useful to consider more general classes of stochastic processes, which are char- 
acterized by more than one parameter, but still preserve some of the good properties of scale 
invariant processes. Therefore, larger classes of stochastic processes have been introduced, for 
example the class of multiplicative harmonizable processes introduced by Borgnat and Flandrin 
[6j, which admits the use of techniques for harmonic analysis. Locally self-similar processes 
(LSSPs) constitute another class of extensions of H-ss processes. This class of processes can 
be used to describe physical systems for which statistical characteristics change slowly in time. 
Based on different situations, several definitions have been raised for LSSPs. Out of the large 
literature, we mention the work of Flandrin et al [T3J, |16j in introduction of LSSPs as Lamperti 
transformation of locally stationary processes in Silverman's sense [39]. Besides, Flandrin and 
Goncalves [18] and Goncalves and Flandrin |19] proposed LSSPs in a sense that the self-similarity 
parameter, H(t), is allowed to vary with time, and discuss their applications. Cavanaugh et al 
[HJ, Courjouly [TTJ, Goncalves and Abry [20] . Kent and Wood [25], Stoev et al [30] and Wang 
et al [45] studied LSSPs in the latter sense, and estimated the local Hurst parameter H(t). In 
addition, Wang et al [35] and Constantine [TO] investigated fitting LSSPs to a geophysical and 
aerothermal turbulence time series respectively. Boufoussi et al [7J studied path properties of 
a class of local asymptotic self-similar process. Moreover, Muniandy and Lim [32] considered 
modeling of LSSPs using multifractional Brownian motion. We also mention the work of Istas 
and Lacauxy [23| on locally self-similar fractional random fields. The road to LSSPs analysis 
was followed especially in time-domain and Hurst parameter estimation. 

In this paper, we consider LSSPs in Silverman's sense [5], [39], and we study their spectral 
analysis. In this sense, a stochastic process is LSSP [5], if its covariance function, Rx{t,s) = 
E(X{t)X* (s)), has the simple form 

R x (t,s) = {ts) H q{\nVts)C x {t/s) := Q(Vte)C x (t/s), 

where Q(t) := t 2H q(lnt), i.e., the covariance function can be decomposed into a covariance 
function of an H-ss process multiplied by a modulatory function. 

Although there is no "universal definition" for nonstationary spectrum, some definitions 
have been proposed for restricted classes. More notable ones are the evolutionary spectrum 
(ES) proposed by Priestley [33] for the class of oscillatory processes, and the Wigner-Villc 
spectrum (WVS) proposed by Martin |30] for the class of harmonizable processes [37] . Because 
of a few important mathematical properties of the WVS over the ES and other nonstationary 
spectrums, such as uniqueness and that the WVS is explicitly defined in terms of the covariance 
function; the WVS has been more desirable than other spectrums [37]. So, we adopt the WVS 
as our definition for the nonstationary spectrum and address the problem of estimating it for 
Gaussian LSSPs. However, the scale invariant property of LSSPs, makes such processes different 
from other nonstationary ones. Therefore, a particular type of the Wigner spectrum should be 
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considered, which is compatible with this property. In this case, for scale invariant signals, a 
bilinear time-dependent extension of the Mellin transform proposed by Marinovic [28] and Altes 
PQ, which will be referred to as the scale invariant Wigner Distribution (SIWD). Flandrin in 
[Li] provided some results of the SIWD and extended it to stochastic processes. 

The scale invariant Wigner spectrum (SIWS) p3] of a process X(t) is defined as an expec- 
tation of the SIWD, and it is a function of time t and frequency £, which reads as 

/oo 
Xit^X^t/^r-^-'dr} (2) 
-oo 

/oo 
-oo 

The usefulness of the SIWS has been emphasized for revealing scale invariant feature in some 
nonstationary processes. By contraction, the SIWS describes the time evolution of Merlin's 
variable in a similar way as the WVS does for the (Fourier) frequency variable [5]. To avoid 
confusion, the classical Wigner- Ville spectrum will be referred to as WVS, and the scale invariant 
spectrum will be referred to as SIWS. 

In this paper, we estimate the SIWS of a Gaussian LSSP, using the mean-square optimal 
kernel method. The optimal kernel method is proposed by Sayeed and Jones [37] for nonsta- 
tionary spectrum estimation; and is based on finding the best estimator that optimizes the 
bias- variance trade-off in the sense of giving minimum mean-square error (MMSE). However, as 
we mentioned, the scale invariant feature of the process and consequently, the spectrum, makes 
the SIWS estimation method different from the other nonstationary spectrum estimations. So, 
we modify the method to be reconcile to our estimation problem. Moreover, as the class of esti- 
mators, we choose Cohen's class counterpart of time-frequency representations (TFR's), which 
is proposed by Flandrin |14] as a generalization of Cohen's class [9] for scale invariant signals. 
Since the Cohen's class counterpart is completely characterized in terms of kernels, the problem 
of estimation is reduced to finding the "best" kernel in the sense of giving the minimum mean- 
square error. The SIWS estimator is obtained in two cases: global case, and local case; where 
in the local case, the optimal kernel is allowed to vary with time and frequency. 

Now, we present an outline of the paper. Next section, concerns the background on self- 
similar, multiplicative harmonizable, circularly symmetric processes. Section 3 treats locally 
self-similar processes, and we study necessary conditions on two functions, that constitute the 
covariance function of a LSSP. In Section 4, we discuss the global and local estimation problems 
of SIWS for Gaussian LSSPs. In Section 5, we introduce two extensions of LSSPs; also, we 
present one example and the performance of the method is studied via simulations. 

2 Preliminaries 

The Fourier transform of a function / £ L 1 (R) is defined and denoted by [41] 

(-77X0 == /(£) == / f(t)e- l2 ^dt. 
A function / : R. 2 -> C is nonnegative definite [1], denoted / G iViVZ)(R 2 ), if 

n 

]T f(t 3 ,t k )z 3 z* k >0 Vfe^cK, {*,•}?=! CC, n>0, (3) 
j,k=i 
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and it is called weakly nonneg ative definite [I], [33], denoted / G WNND(R 2 ), if / is bounded, 
measurable, and 

/ / f(t,s)&(t)d*(s)dtds>0 W£C C (R). (4) 
J Jtr 2 

where C C (IR) denotes the compactly supported continuous functions. If / is continuous and 
bounded, then / G NND(R 2 ) if and only if / G WNND(R 2 ) [3]. For generalized functions, we 
mean by / G NND(S(R 2 )) that / G S'(IR 2 ) and 

(/,tf®#*)>0 VtfGS(R), (5) 

where 5(1^) is the Schwartz space of smooth functions such that a derivative of any order 
multiplied with any polynomial is uniformly bounded, and 5'(lR d ) is its dual, the tempered 
distributions [22], [33]; and the bracket (/, g) = j Rd f(x)g*(x)dx denotes the inner product on 
L 2 (R d ). A function / belongs to NND(R 2 ) if and only if it is the covariance function of a 
mean-square continuous stochastic process |27j . 

A function / : 1R — > C of one real variable is nonnegative definite, denoted / G NND(Wi), if 

n 

f(tj - tk)*jZ*k > V{tj}]=i C R, {zj}?=i CC, n > 0. 

We will use some results about measure theory |36] , The domain of a complex- valued measure 
fi on 1R is a class of subsets of H, which is closed under countable unions and complement, and 
it includes the empty set. We use the Borel cr-algebra, denoted B(1R), as a domain for measures. 
A measure fulfills /i(0) = 0. It is countably additive in the sense of ^(UfcLi ^fe) = Sfe^=i M^fc) 
when {^4fc}fcLi are P a h wise disjoint. The total variation of a measure [i is defined as: 



k=i 



where the supremum is taken over all {A^}^^ C £> such that UfcLi ^fc = ^ an d {^felfcLi are 
pairwise disjoint. The total variation, |//|, is a finite measure if /i is a finite measure. 



2.1 Self-Similar processes 

A stochastic process {X(t),t > 0} is said to be (statistically) self-similar of index H [5] (or 
scale- invariant, denoted 'H-ss'), if for any A > 0, Eq (JJ) is satisfied; i.e., a self-similar process is 
invariant under any renormalized dilation operator Dh,x by a scale factor A [5], |16j . 

A random process {X(t),t > 0} is called wide sense self-similar with parameter H [3BJ, if it 
satisfies the following conditions: 

1. E(X(t)) = \~ H E(X(\t)) for all t, A > 

2. E(X 2 (t)) < oo for each t > 

3. E(X(tx)X*(t 2 )) = \^ 2H E(X(Xt 1 )X*(\t 2 )) for all h,t 2 ,X> 0. 
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Proposition 1 Any wide sense H-ss processes {X(t),t > 0} has a covariance function of the 
form J2J/ 

R x (t,s) = (ts) H C x (-), (6) 
s 

for any t,s> 0, and C x {.) G NND(R + ). 
Proof. See Flandrin [T6l. 



In the case of self-similar processes, the Mellin transform plays, with respect to scaling, a 
role similar to that played by the Fourier transform with respect to shifting [16]. The Mellin 
transform [5], [13] of a function g(t) is defined by 



g(s) := (Mg)(s) = / g(t)t— l dt, (7) 
J o 

where s = /3 + iw, ftweE. 

For a function g of several variables we denote partial Mellin transform with respect to 
variables indexed by j, k, by Adj^Q- 

2.2 Harmonizable, Multiplicative Harmonizable, and Circularly Symmetric 
Gaussian processes 

We present here some definitions about multiplicative harmonizable processes from Flandrin et 
al [5], [6], |16j . A process {X(t),t > 0} is called multiplicative harmonizable if its covariance 
function has the representation 

/oo poo 
/ t H+i27Tl5 s H - i27Ta m(df3,do-), (8) 
-oo J —oo 

where m is a measure on R 2 of bounded total variation, and is called the Mellin spectral 
distribution function. Furthermore, this spectral distribution function satisfies 

poo poo 

m(p,a)= / / t- H - i2 ^- 1 s- H+i2na - 1 R x ( y t,s)dtds. 
Jo Jo 

A necessary and sufficient condition for ([8]) to hold is that \m(d/3, do~)\ < oo, which is 

adopted from Loeve's condition for harmonizability [27] • A multiplicative harmonizable process 
X admits a spectral representation on a Mellin basis as: 

poo ^ 

X(t)= / t H+i2n ^dX(f3), (9) 



oo 



where X is an L 2 (P)-valued measure, called the spectral measure or spectral process. The 
connection between X and m is: 

E(dX((3)dX*(a)) = m(df3,da). 

For multiplicative harmonizable H-ss processes, m has support on the diagonal, and it is a 
nonnegative bounded measure. Therefore, there exists a non-negative bounded measure C such 
that 

poo ^ 

C x (t/s) = / (t/s) i2 ^C(d(3), (10) 
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,C(d/3) 0=<T 

since, by Q and m(d/3, da) = < 

/oo poo ^ ^ 

/ t H+i2 ^s H - i2na E(dX(P)dX*(a)) 
-oo J —oo 

/oo ____ 
(t/ S ) i2 ^C(d/3), 
-oo 

by comparison to ([6]), Eq (fTO]) is achieved. 

Definition 2 For a circularly symmetric or proper process I33f . I38\j X(t), the processes 

{e W X(t)} em27v) 

are identically distributed for all 9 £ [0, 27r). 

According to Grettenberg 's theorem 135)/ . for a circularly symmetric process, E(X(t)) = 
and E(X(t)X{s)) = V(t, s) G R 2 /22F- 

Let Xi, • • • , X n be complex- valued zero-mean jointly Gaussian stochastic variables. Then, 
according to Wick's theorem |24j . 

E(X 1 X n ) = J2U E(X lk X Jk ), (11) 

where the sum is over all partitions of {1, • • • ,n} into disjoint pairs {ik,jk}- Thus for n = 4 we 
have 

S(XiX 2 X 3 X 4 ) = £(XiJSf 2 )£(X 3 X 4 ) + J B(X 1 X 3 ) J B(X 2 X 4 ) + £(XiX 4 )£(X 2 X 3 ). (12) 

The formula (I12p sometimes called Isserlis' theorem is valid for any zero mean complex- valued 

Gaussian stochastic variables Xi,X 2 ,X^,X^ |21j . 

The following fourth-order moment function will be needed later on. 

Remark 3 Let 

fc(ti,ri,t 3 ,Ta) := ^Wiv^^Wv^^fcv^WVv^}- (13) 
Then, from Mty) . for real-valued Gaussian processes, we have that 

k(ti,Ti,t 2 ,T 2 ) = Rx{ti*^,ti/*^)Rx{t2^/^,t 2 /^/v2) 

+R x (hVn, t 2 Vn)Rx{h/Vn, Wv^) 

+ Rx(tiVn,t 2 /V^)Rx(h/Vn^2V^)^ (14) 

and for circularly symmetric Gaussian processes, 

k(h,Ti,t 2 ,T2) = Rx{tiVn,h/ ^)R\(t 2 ^,t 2 / ^) 

+ Rx(hVn,t 2 V^)Rx{h/Vn,t 2 /^ 2 -). (15) 



6 



3 Locally Self-Similar Processes 

Locally self-similar processes (LSSPs) in Silverman's sense [39], were introduced by Flandrin 
|16] as an extension of H-ss processes. In this Section, we study properties of LSSPs in time and 
time-frequency domain. 

Definition 4 A locally self-similar process (LSSP) J2J/ ; 11 6\/ . is a complex-valued stochastic pro- 
cess {X(t),t > 0} whose covariance function Rx(t,s) = E{X{t)X* (s)), has the form 

R x (t, s) = (ts) H q(ln VTs)C x (-) (16) 

s 

where q and Cx ere complex-valued functions, q must have a constant sign which we assume pos- 
itive, and Cx G NND(1R, + ). If we define Q(t) := t 2H q(lnt), and the coordinate transformation 
on R +2 as 

kit, s) := (t^/s, t/y/s) <^ fc _1 (*, s) = {Vts~, -), (17) 

s 

then {Tty) can be written as 

Rx(t,s) = Q(Vts~)C x (-), (18) 

s 

we may write this relation as 

R x {t,s) = Q®C x °k- 1 {t,s). 

IfRx G NND(R+ 2 ) and R X = Q®C X ° kT x , we write (Q,C X ) G L55(R+). 

We will always understand that the stochastic process is nonzero, i.e., Q, Cx ^ 0. Let (Q, Cx) £ 
L5'5(]R + ), and let X(t) be a locally self-similar process with covariance function Q <g) Cx ° k~ l . 
Then 

Q(t)C x {l) =E\X(t)\ 2 
describes the variation of the process over time, and 

Cx(r)/C x (l) = E{X{t^)X*{t/^))/E\X{t)\\ 

describes the "local covariance" (at time t, for any t G R+ such that E\X(t)\ 2 / ), of the 
stochastic process X(t). Since X is assumed to be nonzero, we have Cx(l) 7^ 0, and we may 
without loss assume that 

C x (l) = l for (Q,C X )€LSS(R + ). (19) 
Thus, Q(t) > for all t G R + . 

We also assume that Q, Cx are continuous. Continuity of Q and Cx is equivalent to mean- 
square continuity of the process X(t). In fact, everywhere continuity of Rx is implied by 
continuity on the diagonal [27]. Thus, if Q is continuous everywhere and Cx is continuous in 1, 
then Rx is continuous everywhere, and the process is mean-square continuous. 

The Cauchy-Schwarz inequality for the Hilbert space L 2 (Q) gives 

\Rx{t,s)\ 2 < R x (t,t)R x {s,s), 
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and is a necessary condition for a function Rx to be the covariance function of a stochastic 
process [27]. If (Q,Cx) £ L55(R + ), we have due to the normalization (fT9j) which implies 
Q > 0, 

Q 2 (t)|Cx(r)| 2 < Q(ty/7)Q(t/y/7), t,re R + . 
By restriction to i = 1, we obtain the inequality Q 2 (1)|Cx(t 2 )| 2 < Q(t)Q(1/t), or equivalently 

Q 2 (1)|C x (t 2 )| 2 <Q(t)Q(1/t), (20) 

which says that Cx is necessarily bounded in terms of Q(t)Q(1/t). 

We are interested in conditions on Q and Cx that are necessary and sufficient for (Q, Cx) £ 
LiS\S(<S(]R + )), which is equivalent to (Q,Cx) £ LS'S'(R + ) when Q an d Cx are continuous 
and polynomially bounded. The following simple lemma expresses the definition (Q,Cx) £ 
LSS(M 1 + ) in terms of time-frequency analysis. 

Lemma 5 IfQ,C x £ S(R + ) then 

poo poo 

(Q®C X ) £ LSS(S(R + )) / tQ{t)C x (€)Wf(t/to,S-So)dtdt>0, f £ S(R+) (21) 

JO J-oo 

where Wf(t,£) = J °° f{ty/r)f*(t/y/r)T~ %2 ' K ^~ 1 dT, is the scale invariant Wigner distribution. 

Proof. By (Q, C x ) £ LSS(«S(IR+)), we mean R x = Q^C^k" 1 G iViVL>(R +2 ) with continuous 
and polynomially bounded functions Q and Cx- Then, by ((5J), (Q<8> C x o /c" 1 , (8) i?*) > 0, Vi? £ 
«S(R+). Thus, 

/*oo /*oo 

(Q®C X ok- 1 ,^®^*) = / Q^CxoAT 1 ^) ^* ®$(u,v)dudv 

Jo Jo 



oo />oo 



J 



Q( v / ^)C x (n/t;)i?*('u)i?(v)d^ 



/•00 f'OO roo 

= / / tQ(t)C x (0{ ^(t^Wt/^y^dTjdtdt^O 

Jo J-oo Jo 

where t := ^/uv,t := u/v and Cx '■= M Cx- By real-valuedness of we have that 

roo roo 

{q®c x o ArV<g>r) = / / tQ(t)Cx(OW0(*,O^t > o. 

Jo J-oo 

Let / G S(R+) and set 0(t) := t i2 ^°f{±) G S(R+), then = W>(t/t ,£ - Co), and (J2TJ| 

is satisfied. 

□ 
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3.1 Multiplicative harmonizable locally self-similar processes 

In the remainder of this section, we study locally self-similar processes that are multiplicative 
harmonizable. Let there exists a bounded measure C such that ()10p holds, and also there exists 
a bounded measure Q such that Q(t) = JZo t2H+i2 ^Q( d O- Th en for a LSSP, © takes the 
form 

/oo poo ^ ^ 

/ (^s) 2H+i2 ^(t/s) i2 ™Q{dOC{duj) 
-oo J —oo 

oo roo 



oo J — oo 



By the change of variables j3 := u + £/2 and a := cu — £/2, we have that 

J —oo J—oo " 

By comparison to ([8]), it follows that the spectral distribution function, m, can be written as a 
product of bounded measures C and Q : 

m (^a) = Q(/3- C 7)C(^). 
4 The Scale Invariant Wigner Spectrum Estimation 

In this section, we study the SIWS estimation of a zero mean real-valued, and complex-valued 
Gaussian LSSP defined on a positive half of a continuous time axis. As the class of estimators, 
we use the Cohen's Class Counterpart |14j of time-frequency representations (TFRs) for scale 
invariant signals. The class of estimators, closely parallels the conventional Cohen's class [9] and 
shares with it some of its most interesting properties, namely those concerning the usefulness and 
versatility of distributions associated to separable smoothing functions [14j . We study estimators 
in this class that are optimal in the sense of giving minimum mean-square error (MMSE). 

4.1 The Scale Invariant Wigner Spectrum 

Given a scale invariant signal x defined on R + , the scale invariant Wigner distribution (SIWD) 
HI, [2S], is defined by 

roo 

W x (t,0= x^x^t/^r-^-'dr (22) 
J o 

= M 2 (x ® x* o fc)(t,r), 

which can be interpreted as a distribution of the signal's energy over the time-frequency domain 
£ 1R + x 1R. However, this interpretation is bothered by the fact that W x is rarely non- 
negative everywhere. As a remedy, one can smooth W x with a kernel according to W x $, 
which is called a time-frequency representation in Cohen's class counterpart determined by P x 

EH 
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If X is a stochastic process, then W x is called the scale invariant Wigner process, and is 
denoted by Wx- In this case, the integral (|22p is a stochastic mean-square Riemann integral 
[27] ; and the scale invariant Wigner spectrum (SIWS) is defined by the expectation of Wx' 

poo 

W E , x (t,0=E{ Xit^X^t/^r-^-'dr} (23) 
J o 

poo 

Jo 

= (M T Rx)(i2n£) 

where Rx is the covariance function of a process {X(t),t > 0} [5], |14j . 

Proposition 6 In the case of multiplicative harmonizable processes, the SIWS can be equiva- 
lently expressed in terms of the spectral covariance function, m(/3,o~), as 

/oo 
t 2H + i2 * e m(Z + e/2,Z-e/2)de. (24) 
-oo 

Proof. See Al. 

The scale invariant ambiguity function (SIAF) [TJ], can be obtained from the SIWD by 
Mellin duality: 

poo 

A x (6,t)= / x(t^)x*(t/^)r i27r0 - 1 dt (25) 



o 



= MiM^W x (t,C), 

and for a stochastic process X, the Expected scale invariant ambiguity function (ESIAF) is 
defined by: 

poo 

A E)X (0,t) := E{A x (9,t)} = / R x (tV^,t/^)r^ e - l dt 

Jo 

= (MtRx)(i27T0) 

= M 1 Mz 1 W E , x (t,0. (26) 



For a locally self-similar process, by (|18l) . 

Rx(tV^,t/^) = Q(t)C x (T), 
so by (j23|) . one can easily verify that the SIWS of a LSSP can be represented by: 

poo 

W E ,x(t,0= Q(t)C x (r)T~^ l dT 
Jo 

= Q{t)(MC x )(i2^), (27) 

and the ESIAF for this class is 

poo 

Jo 

= C x (T)(MQ)(i2ir9). (28) 
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4.2 The Cohen's Class Counterpart 



The Cohen's class counterpart of bilinear time-frequency representations (TFRs) of scale invari- 
ant signals is represented as 

poo poo ± J 

Px(t,0= / Wsi-^-riWs^—dri, (29) 

J-oo JO s s 

where <]? is the 2-D kernel that completely characterized the particular TFR P x , and W x is the 
SIWD [J2]. Eq ([29]) can equivalently be represented as [T4"] 

/oo poo 
/ A^e^ie^t^r-^^dTde (30) 
-oo JO 

= M^M 2 {A x {e,T)<t>{e,T)}, 

where </>(0,t) = f™^ J™ ' £)i _i27re_ V 27r Qid£. This class, closely parallels the Cohen's class 
[9] , and it shares with it some of its most interesting properties [H] . 

If x denotes a realization of a scale invariant process {X(t),t > 0}, then the integrals defined 
in (|29p and (|30p become stochastic integrals and will be interpreted as mean-square (m.s.) 
integrals [27], [37]. We assume that the kernel (ft is chosen such that the existence of the SIWD 
as a m.s. integral, implies the existence of Px as a m.s. integral. 

4.3 The SIWS Estimation 

Here, we study optimal kernel estimation of the SIWS, for Gaussian LSSPs in two cases: (i) the 
local case in which the kernel (p is allowed to vary with time and frequency in order to better 
track the nonstationary structure of the process, and (ii) the global case, where the kernel <j) is 
assumed to be independent of time and frequency. 



4.3.1 The Local Estimation 

Let {X(t),t > 0} be a zero mean real-valued, or complex-valued circularly symmetric, Gaussian 
process. Then, (f22l) . (f25l) and ([30]) will be the stochastic Riemann integrals, denoted by Wx, 
Ax and Px- The stochastic integrals exist for all argument values [27], and are second order 
stochastic processes [43]. Px(t,£,), by definition, is a member of Cohen's class counterpart of 
TFRs. 

Now, we consider the local SIWS estimation problem. We use MMSE estimator where the 
optimality criterion of the estimation will be the minimization of the mean-square value of the 
estimation error. So, our objective is finding the optimal kernel <f>^fP, which minimizes the 
mean-square error 

J(4>^) = E{\P x (t,0 - W E , x (t,0\ 2 } (31) 

for each (t, £). This is a problem of linear MMSE estimation. As we mentioned, in the local 
case, to better track the nonstationary structure of the process, the kernel may depend on time 
and frequency, and the superscript (£,£) denotes this possible dependence. 
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First, note that for each value of (i,£)> {Px(*,£);^ £ L 2 (R x R + )} belongs to a Hilbert space 
% of second-order random variables. The scale invariant ambiguity function Ax(9,t) generates 
a subspace Sx of % as defined by Px(i,£) in (f30j) . ref [37] : 

5x = { J Px(t,e):<AG^ 2 (HxR+)}. (32) 

In the local problem (f3Tj) . the orthogonal projection of W^xC^O onto Sx is desired. So, using 
the expression for Px(t,£) in (I30j) and by orthogonality principle, Px(t,£) achieves minimum 
mean-square error if and only if [37] 

E{[P x (t,O-W Et x(t,O]Ax(0',r')}=0, for all (9' ,r'), (33) 

where A* x is a complex conjugate of Ax- 
Expressing Px as in (f30|) . we can write (|33|) as 



OO fCO 



E{Ax{e,T)A\{e'y)}^\e,T)^ e T- l2 <- 1 dTde 

oo JO 

= ^,x(U)M4(^V')} 1 for all (9', r'). (34) 

The above linear equation characterizes the locally optimal kernel The linear equation is 

of the form 

A [W 4> = W EyX (t,0b, (35) 

where A : -L 2 (R x R+) ->■ L 2 (R+ x R) is a linear operator, and 6 = EA* X G L 2 (R x R+). Thus, 
the solution can be written as 

^=W E ,x^OA\ t>0 b, (36) 
where the superscript f denotes the pseudo-inverse. By insertion of f)35p into (l30|) . we have that 

/oo roc 
/ Ax(0,r)05?(^T)t i27r0 T- i2 ^- 1 drd0 
- oo JO 

/OO /'OO 
/ A x (9, T){A\ t M™T- i2 ^drdO 
-oo JO 

= W B! x(t,e)^(t,C), (37) 
this means that, for each value of (i, £), the optimal estimate of W Ej x(t,0 from a realization 
of a scale invariant process X is W^x^jO itself, multiplied by a TFR of x generated with the 
kernel AL ^b. 

4.3.2 The global estimation 

Now, if we consider (p op t to be independent of (£,£), then by multiplying both sides of (j34|) with 
t _; 2 7r0'-i( T /y27r^ and i nte grating over (£, £), we obtain the global optimal ambiguity-domain 
kernel as: 

^ t)=1 40w^ Iu ^ t) ' (38) 

where Ijj(9,t) denotes the indicator function for the open set 

U = {{9,t);E\A x {9,t)\ 2 > 0} C R x R+. 
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Proof. SeeA2. 

The optimal time-frequency kernel is computed by 

<5>opt(t,Z) =M^M 2 M0,t). (39) 

For circularly symmetric Gaussian processes, using (|15p and (j26|) we obtain 

E\A x (6,t)\ 2 = \A e , x (6,t)\ 2 + D^e^), (40) 

where 

fOQ f'OO 

£>i(0,t)=/ / Rx(tiV^,t 2 ^)R* x (ti/V^,t2/V^)ti t2ne - l tf^- 1 dt 1 dt 2 . (41) 
io JO 

In the case of real- valued processes, using {H}, we get 

S|ylx(^r)| 2 = \A E , X (6,T)\ 2 + D 1 (e,r) + D 2 (0,t) (42) 

where 

f'OO f'OO 

D 2 (9,t)= / / R x (t 1 V^,t 2 /^)Rx(ti/V^,t 2 V^K l2ne - l tf 7re - 1 dt 1 dt 2 . (43) 
For validity of (gDJ) and (g2|), see A3. 



Proposition 7 The optimal kernel for a circularly symmetric Gaussian LSSP, is 

6 (f)T) = \C x (r)\ 2 \(MQ)(i2^)\ 2 

9aptK ' ' |Cx(r)| 2 |(^g)(i2vr^)| 2 + (^|C7x| 2 )(i27r0) (Q®Q*)(r)' 1 ' 

where (Q ® Q*)(t) := J °° Q{uy/r)Q* (u/ \/r)^ ■ To obtain the optimal kernel for a real-valued 
Gaussian LSSP, we have to take into account the third term D 2 (8,t) into the denominator. It 
is 

f°° dii 

D 2 (e,r) = A Cx (0,r 2 ) Q 2 (u)-, (45) 
Jo u 

where, A Cx (9,t 2 ) := J™ C x {vt)C x (v /^v'^^dv. 
Proof. SeeA4. 

5 The Optimal Kernel for Some Generalizations of LSSPs 

In order to generalize the LSSP definition, we introduce two classes of locally self-similar chrip 
process, and Multicomponent Locally Self-Similar process. 
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Fi gure 1" Example of covariance matrices R,x of LSSPs with H = 0.5 and different parameter values. Top left: LSSP: 
c = 1.1. Top right: LSSP, c = 30. Bottom left: MLSSP, ci = 1.1, c 2 = 30. Bottom right: LSSCP, c = 1.1, a = 2, b = -2. 

5.1 Locally Self-Similar Circularly Symmetric Gaussian Chrip Process 

Definition 8 A locally self-similar chrip process (LSSCP) is a Gaussian circularly symmetric 
process with covariance function of the form 

R z (t, s) = Q(Vte)C x (t/s)l a>b (t, s), (46) 

where l a ^ is defined by 

la,b(t, s) = (t/s) ia ^ nVrs - b \ a, b G R. (47) 

The constant a determines the chrip frequency and b the start of the chrip frequency. 

If X(t) is a LSSP with covariance Q®C X ° AT 1 , then Z(t) := X(t)t ia( - ln ^t-b) i. e . ; z equals X 
times a chrip, has covariance Rz- We have Rz (ty/r,t/ y/r) = Q{t)Cx{T~)T ia ^ nt ~ b \ which shows 
that C(r) is modulated by the chrip T m ( lnt ~ b ) . 

If one multiplies a signal with a chrip according to f(t) := fo(t)t l ( a / 2 ^ lnt ,a £ R, then the 
SIWD changes to W{t,£) -> W(t, f - £ hit): 

/>oo 

W>(t,0 = / fit^fit/^T^^dr 
Jo 
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Figure 2: Contour plots of covariance matrices Rx of LSSPs for H = 0.5 and different parameter values (corresponding 
to Figure 1). Top left: LSSP, c = 1.1. Top right: LSSP, c = 30. Bottom left: MLLSP, c\ = 1.1, c 2 = 30. Bottom right: 
LSSCP, c = 1.1, a = 2, b = -2. 



roc 

Jo 

= Wf (t,t-£lnt). 

The next theorem states that the optimal kernel for an LSSCP is obtained by a transforma- 
tion of the kernel of the corresponding LSSP. 

Theorem 9 Suppose Cx and Q define a circular symmetric Gaussian LSSP with covariance 
given by $18\). Then, the optimal time-frequency kernel for a LSSCP process with chrip param- 
eters a and b, denoted <& c .opt, is 

<S>c.opt{t, = <S>o P t{t, £ - ^ In t) (48) 

where Q op t is the optimal kernel of the corresponding LSSP (i.e., a = b = 0), defined by /139\) 
and (44\ ). 

Proof. SeeA5. 
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Figure 3: Contour plot of covariance matrices Rx of LSSPs for different parameter values c and H. Left column: LSSP, 
c = 1.1, and H ranging from (top) H = 0.05 to (bottom) H = 0.95. Middle column: LSSP, c = 5, and H varies from 
H = 0.05 to H = 0.95. Right column: LSSP, c = 30. Each column shows that the covariance shape for a constant c is 
identical, but as H increases, the slope of Rx decreases. Also, each row shows that, with the same Hurst parameter H, as 
the parameter c increases, the peaks in sides are inclined into the peak in the middle. 

5.2 Multicomponent Locally Self-Similar Circularly Symmetric Gaussian Pro- 
cess 

Definition 10 A multicomponent LSSP (MLSSP) is a process whose covariance has the form 

oo 

Rx{t,s) = Y J QA^)Cx 3 {t/s), (49) 
j'=i 

where each term Qj <H> Cxj ° k~ l is the covariance function of a LSSP. 
By insertion of gS} into (123]) and M, the SIWS and ESIAF are 



52Qj(t)(MC Xj )(i2rt), 



3=1 



A 



E,X{V,T) 



POO °° 

Jo , =1 

CO 

3=1 
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c=1.1,a=2,b=-2 





I _ 2 



Figure 4: Sample realizations of LSSPs with H = 0.5 corresponding to the covariance function Rxit, s) in Figure 1. Top 
left: LSSP, c = 1.1. Top right: LSSP, c = 5. Bottom left: LSSCP, c = 1.1, a = 2, b = -2 (real-valued part only), and 
Bottom right: MLSSP, ci = 1.1, C2 = 30, which is a combination of LSSPs with c = 1.1 and c = 30. 

Restricting to the circularly symmetric case; and using (j40j) and (|4ip . the optimal ambiguity 
domain kernel for circularly symmetric MLSSP is derived as 



ET=iC Xj (TKMQ 3 )(i2Ke)\< 



Ef=i C Xj {T){MQ 5 )(i2ire)\* + Y^ k=l M{C Xj C Xh )^e) (Q j © Ql)(r) ' 
where the second term of the denominator is 



(50) 



di(o,t) = J™ J™ f; Q j (vW)Q^(^)c-x,(|)c'^(|)tr^- 1 4 2 ^ 1 ^i^2 

j,k=l 

roc 7 

As a generalization of LSSCPs and MLSSPs, we can introduce a multicomponent locally 
self-similar chrip process (MLSSCP) whose covariance function has the form 



R X ( t , s ) = Y,Qj(^)Cx 3 (t/s)(t/ s y a > {lnVFs - b > ) , 



(51) 



where each term are LSSCP covariances with individual constants dj and by A MLSSCP has a 
covariance that is a sum of covariances with chrip behavior of various localization in time and 
chrip constants. 

By insertion of (f5"T|) into (f26j) and (|4T|) . we have 
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* opt (6 |t ], c = 1.1,H=0.1 



+ opt (8, t), c = 1.1, H=0.5 



♦ „ pt (e,i), c = 1.1,H=0.9 




Figure 5: The mean-square error optimal ambiguity domain kernel corresponding to the LSSP with covariance function 
in Figure 1, with c = 1.1 Left: 0^f t SP (0, -r), H = 0.1. Middle: tfffi ,F '(0 , r) , H = 0.5. Right: <f>^ t s p {9 , r) , H = 0.9. 



A E ,x(e,T) 



f'OO °° 

00 roc 
*£C Xj (T)T- ia ^ / Qj(t)t 

,=i h 



-i2ir0+ia; Inr— 1 



poo poo ^ 

Di(0,r)= / / E Q j (V^)QU\l—)c Xj mc Xk Cr) 

Jo Jo f^ x V r t 2 i 2 



x ( — ) ia J ( ln V*i*2T-bj)-iafc (In Ht 2 /T-b k ) £-12^0-1^2^0—1^^ 
*2 



E 



i,fe=i 



o Jo 



-i27r{0-|i(ln U V?-6,)+|i(ln^-6 fc )}-l *i 



civ. 



Thus, we have that 



OO n QQ j 

£>i(0,r) = V(M|Cx/)(^-^lnr)) / QjiWr)^/^ — 
~^ ^ jo u 



+ E 



00 ciu 
Qj (uVt)QI (u/y/r)g(6, u) — , 

u 



(52) 



(53) 
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where 

g(9,u) = M(C X] (v)C* Xk ){9 - g(ln V^- bj) + g(ln-^ - 6*)). 



By insertion of (I52p and ([53]) into (138p . the optimal ambiguity domain kernel for a circularly 
symmetric Gaussian MLSSCP is obtained. 

5.3 Discrete-time Examples 

Here we study an example of a circularly symmetric Gaussian LSSP, LSSCP and MLSSP with 
covariance functions (fl8l) . (j46l) and (f49l) respectively. We consider 

Q(r)=T 2H -^ nT , C x (t) = r-f lnr ,c> 1. (54) 

For a given Q and Cx, the covariance matrices, Rx, of the processes are computed, and depicted 
in Figure 1, for H = 0.5, c = 1.1, c = 30; and the corresponding contour plots are shown in 
Figure 2. Also, the contour plots of the LSSP for three different parameter values of H ranging 
from H = 0.05 to H = 0.95, and c = 1.1, c = 5 and c = 30 are displayed in Figure 3. 

The sample paths of the circularly symmetric Gaussian LSSP, LSSCP and MLSSP, corre- 
sponding to the covariance matrices of Figure 1, are simulated and displayed in Figure 4 with 
H = 0.5. The realizations are simulated from the covariance matrices Rx, according to 

X = Lu, (55) 

where u is a realization of a real- values white Gaussian zero mean stochastic process with variance 
one, and the matrix L is related to the covariance matrix Rx as 

R x (t, s) = E(X(t)X*(s)) = L E(u(t)u*{s)) L* = LL*. (56) 

The optimal ambiguity domain kernel for a circularly symmetric Gaussian LSSP, LSSCP 
and MLSSP are obtained in (f4"4"j) . (f63|) and ([50]) . where for given Q and Cx in (f54"|) are computed 
as 



i MLSSP i 



ALSSP (f) s = 1 ( , 7) 

9 ° pt K ' ' 1 + c-1/2 e (l-l/c)(2ner+AHi(2^) r (^±) Mr ' 

b LSSCP,Q n = 1 / 58 n 

° P * K ' ' l _j_ c -l/2 e (l-l/c)(27r6»-alnr) 2 +4Hj(27re-alnr) T ( £ T i ) lnr ' 

V /^ e -^( 2 ^) 2 +(^'+^) 2 r (H 3 -H k )-|ln, 

V(l + EWcj+C :! ^ ^ )], (59) 



" l E . c i(2H,-^)V-^Inr|2 



^. r~2~ c -^ r (2^alnr) 2 +(H J +H fc ) 2 ^-^J-l l nT 



<t%t SSCP {P,T)= (l/(l + ^' fc V^ ., (6Q) 

P l'\ |^. e i(2^-i(2 7 re-alnr))2 T -ilnr|2 ' 

Eq (|60p is the optimal kernel of a MLSSCP which is derived as a combination of the optimal 
kernels for LSSCP and MLSSP. For validity of (|57p - (l60p . see A6. The optimal ambiguity domain 
kernels of the LSSP, Eq (|57p . for c = 1.1 and three different Hurst parameters ranging from 
H = 0.1 to H = 0.9 are plotted in Figure 5. 
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Time Time 



Figure 6: Left: True SIWS of the Gaussian LSSP. Right: The SIWD of the Gaussian LSSP. 

5.4 Simulation 

Now, the predominance of the proposed method for the SIWS estimation of a Gaussian LSSP 
is investigated over the classical optimal kernel WVS estimation. To this end, the performance 
of the globally optimal kernel estimator is compared with that of the classical WVS. 

Let X be a Gaussian LSSP with the covariance matrix (|18h where Q and Cx are identified 
by ((SU). We assume that H = 0.5 and c = 1.1, so 

R x (t, s ) = ^-^^{t/syi *(*/->. 

The SIWS of a LSSP is given by ([27]) . and is computed for given Q and Cx , Figure 6 (Left). The 
SIWS is estimated from the realization X, using relation (|30p . where the optimal kernel 4>{0,t) 
is computed in (|57p . Figure 7 (Left). Clearly, the optimal kernel SIWS estimation gives a much 
more accurate estimate for the true SIWS, than the classical optimal kernel WVS estimation 
Figure 7 (Right). 

The simulation results show that, although LSSPs are a subclass of nonstationary processes, 
but the scale invariant property exists in such processes, makes them different from the other 
nonstationary ones; and ordinary nonstationary spectrums and estimation methods may not be 
applicable for TF analysis of such processes. So, some special tools should be considered, which 
are compatible with the scale invariant property. 
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Fi gure 1\ Left: The optimal kernel SIWS estimation of a Gaussian LSSP. Right: The classical optimal WVS estimation. 

6 Conclusions 

In time- frequency analysis of locally self-similar processes (LSSPs), we obtained the scale in- 
variant Wigner spectrum (SIWS) estimation by using the Cohen's class counterpart of time- 
frequency representations, which is compatible with the scale invariant property. By this, we 
provided a better estimation than the classical Wigner- Ville spectrum. By introducing the op- 
timal kernel for the SIWS, which minimizes the mean-square error, our estimation has been 
modified and provided a close estimation to the true SIWS. 

Although we restrict ourselves to Gaussian LSSPs, but the MMSE optimal kernel can be 
computed for the SIWS estimation of non-Gaussian LSSPs. In such problems, the optimal kernel 
would be as a linear combination of the optimal kernels of Gaussian LSSPs. 

7 Appendix 

Al: Proof of proposition 6: If X{t) is a multiplicative harmonizable process, then from (JSj) we 
have that: 

/OO /'OO 
/ (t^f +i2nP (t/^) H - l2na m(P, a)d(3da. 
-oo </ — oo 

By insertion of the above relation into the definition of the SIWS , Eq (|23[) . 

poo poo poo f , 

W E , x (t,0= / / t 2H+ ^^T^^T- i2 <- l m{p,a)dpdadT 

Jo J — oo J — oo 
poo poo poo 

= / t 2H+i2vd T- i2 ^-^- l m{f + 9 /2,f - 6 /2)d6dfdr, 

J J — oo J —oo 

where 6 := j3 — a and / := ^rp. Using relation J °° T~ i2% ^~'^~^dT = 5(£ — /), where S(.) is a 
dirac delta function, (|24|) is achieved. 
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□ 



A2: Validity of Eq (|5gj) : By multiplying both sides of ([Ml) with i-i27r0'-i( T /)i27r£ and integrating 
over (t, £), the l.h.s leads to: 

/OO /"OO /'OO /"OO 
/ / / E{^ x (0,r)^(0',r / )}^ f (0,r)t-^( e '- e )- 1 (^)- i2 ^-^^ 
-oo Jo J-oo io r r 



= / / E{A x {e,T)A* x {9'y)}<t> opt {e,T)8{d' -e)5{hxr-\nr')— d6 

J-oo JO T 

= E\A x (6,T)\ 2 opt (6,T), 

where / °° t^^dt = 5{u) and v~ i2 ^d^ = S(hxv). By insertion of W E ,x = M^MzAejc 
into the r.h.s. of (IMP, we have that: 



/oo /■oo poo poo j 
/ / / A E ^6,r)E{AWy)}t- t2 ^- e) -\^r 2< -d9dtd^ = \A E , x (9,r)\ 2 . 
-oo JO J-oo JO T T 

□ 

A3: Validity of (00} and (02}: By (p5j). 

^|^(e,r)| 2 = S|A x (0,r)^(0,r)| 

( /"OO /"OO ^1 

= ^{y / xih^x^h/^x^h^x^/^t^-hf^dhd^i 

roo roc 

= / / EiXih^X^h/^X^h^X^/^t^-Hf^dhdh. 

JO JO 

The integral inside the expectation operator is a stochastic integral. The interchange of expec- 
tation and integration is justified if the above-mentioned stochastic integral exists in m.s. sense 
[37]. For a circularly symmetric Gaussian process, using (|15h . we have that 



poo poo 

E\A x (9,r)\ 2 = / Rx(h^,t 1 /^)Rx(t2V^,t2/V^)ti i2ne - 1 4 27re - 1 dt 1 dt 2 
Jo Jo 

rOC rOC 

+ / / Rx(tiV^,t 2 V^)R x (h/V^,t2/V^)ti l2n0 - 1 t¥ n0 - 1 dt 1 dt 2 . (61) 

JO JO 

By the ESIAF, ([26]) . the first integral is \A e ^ x (9,t)\ 2 , and we define the latter one as D\(9,t). 
For the case of real- valued processes, using ([14]) . we have an extra term ([43]) in E\A X {9 ,t)\ 2 . 

□ 

A4: Proof of Proposition 7: For a circularly symmetric Gaussian processes, the denominator of 
(jMD is obtained from gD}. By (ggj), A EjX (9,t) = C x (T)(MQ)(i2ir9); and £>i(0,t) is obtained 
from gU by R x (t,s) = Q(\fts)C x (t/ s) as 
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D x (0,t) = Q{V^)Cx{h/t 2 )Q%^^)C x (h/t 2 )t^ e - 1 tf^dhdh, 



let u := \Jt\ti and v := t\/t 2 - Then, 



D 1 (6,t)= I I Q{u^)Q*(u/^)C x (v)C x {v){u^r^\u/^Y 2 ^ l -dudv 
Jo Jo v 

Q(uV^)Q* {u/^)C x {v)C* x {v)v- i2 ^- l dv^ 




JO 

oo 



\C x (v)\'v 



2„ .-127T0-1 







/ Q{uy/?)Q* {u/y/P) — 

Jo u 



f°° dii 

(M\C x \ 2 )(i27T9) / Q(«V^)Q*(«/V^)— • (62) 
Jo n 



For real- valued Gaussian LSSPs, by (jl3 



oo /-oo 



D 2 (0,t)= I I Q 2 (V^)Cx( t ^T)Cx(^)q l2 ^ 1 tf^ 1 dt 1 dt 2 , 

Jo Jo *2 *2T" 

oo /-oo J 

Q 2 fttjCTx (wr)Cx (v/t^-* 2 * ' 1 — d« 
o Jo « 

Q 2 {u)— / CxHCx^/t).-* 9 - 1 ^, 
o u Jo 



and (H5l) is obtained. 



A5: Proof of Theorem 9: By insertion of (|46l) into (j26H and ()4ip . we have that 

/■oo 

A BjX (0,t) = / Q(t)C x (r)T ia( ' lnt - b h- i27re - 1 dt 
Jo 

poo 

= T~ iab C x (T) / Q(t)t" i2n ^"^ lnTS >~ 1 dt 







T- iab C x (T)(MQ)(i2^{9 - — lnr)), 

Z7T 



□ 



JO JO *2 «2 V T t 2 E2 



JO 



Q(«^Q*(«/^Cx(t;)C^(t;)« < ^^- 6 )«- ia(1,l ^- 6) t;- i2 ^- 1 — d«, 

it 
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where u := y/t\ti and v := f*. So, 

f'OO f'OO 

D 1 (9,t)= iCxiv^v-W-ft^^dv x / Q{u^)Q*(u/^) 
Jo Jo 

= (M\C x \ 2 )(i2ir(6 - ^lnr)) f°° Q( U ^)Q*( U /^) — 
2ir J u 



u 



(M\C x \ 2 )(i2ir(8 - ^-lnr)) (Q®Q*)(r). 



Thus, using ([38]) , (|40p and (|44p . the optimal ambiguity domain kernel for the circularly symmetric 
Gaussian LSSCP, denoted (frc.opt, is given by: 

, , ) = \C x (r)\ 2 \(MQ)(i2n(e-^lnr))\ 2 

<t>C.O P t[ ,V | C - x(T) |2| ( ^Q )(i2 ^-|.l nT ))|2 + ( > (| ax |2 )(i27r(0 _^ lnT))( Q @ Q* )(r) - 

(63) 

By ([11]), it is clear that <j> c .opt(Q \ T ) = 4>opt(6 — ~§^^ YlT i T ) where (j) op t is the optimal kernel of 
the LSSP corresponding to a = b = 0. Now it follows from = M^_ l Ai2 ( t ) that 

poo /*oo 

<S>c.opt(t,0= / <t>c.opt(e,T)? 2 " e T- i2 <- l dOdT 

Jo J-oo 

poo poo 

= / / ^-^lnT,^ 2 ^- 42 ^- 1 ^, 

JO J-oo 27T 

let 0' :=6- g- lnr, then 



OO r-oo 



Jo J-oo 



□ 



A6: Validity of ([5T ]) - ([6T)]) : By insertion of functions Q and Cx defined in (JMD in Eq pi]) , the 
optimal kernel (|57p is achieved. 



/'OO /*oo 

(MQ)(i2TT0)= Q{t)t- i2n0 - 1 dt= t^-^-^^dt 
Jo Jo 

[•OO 

JO 

where u := Int. So, the numerator of (|44p is obtained as: 

Us,x(^,r)| 2 = \C x (T)\ 2 \(MQ)(i2TTd)\ 2 = 2vr T ^(^)^ e ( 2 ^ i2rf ) 2 , (64) 



Also, 



r>00 

(M\C x \ 2 )(i2ir9) = I T~^ nT T~ i2 ^ e ~ l dr 
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OO POO 



§lnr-i27r6»-l. 



T * 



and 



" Q(uyfiF)Q*{u/yfiF)— = r \uV^) 2H - 1 > ln ^H^) 2H - 1 > HU/ ^ ) - 

Jo u Jo V T u 



roo 



Thus, 

/>oo 

D 1 (9,t) = (M\C x \ 2 )(i27r6) / Q(u>/7)Q*(u/V^) 

JO 



it y^c 

By insertion of (|64p and ()66f) in (|44j) . we have that 

2?r r -(c/4)lnr e (2H-i27r9) 2 



27F e -|(7r9) 2 +4^ r _i lnr _ 



27r r ~(c/4)lnr e (2/f-i27re)2 , 2vr g -| (^0)2+4^2 i lnr 



1 + c -l/2 e (l-l/c)(27re)2+87ri0/f r ( £ T L )lnT 

Validity of Eg (p)): 



/>oo />oo 
{MQ)(i2TT(0-— lnr)) = / Q(t)t~ i2 ^ -^ ln ^" 1 dt = / ^-ilnt-^-^ Inr)-!^ 

= / e -|« 2 +«(2^i2^^1nr)) du = ^ e (2H-i27r(0-^lnr))2/2_ 

Jo 

Then, using (j67H and |Cx(t)| 2 = t~( c / 4 ) lnT , the numerator of (|63H becomes 
U £X (^ T )| 2 = |C x (r)| 2 |(^Q)(i2^-^lnr))| 2 = 27r T -(c/4)inr e (2H-i^-^inr))« 

' 27T 

For the denominator, we have that 

(M\C x \ 2 )(i2ir(9-—liiT))= / t -i^t t -i2^0-alnr)-l dt 
27r Jo 

Jo 

=jT 

./o 

So, by ([65]) and (HD, 



e 4 

'o 



s(u 2 + ^2,(«-i!nr)«) rfu 



Di(0,r) = (.M|C x | 2 )(i2^(# - ^- In r))(Q©Q*)(r) = 2^ C - 1 / 2 e 4H2 - 4 ^-^ ln 

27T 
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By insertion of ([68]) an d (|70p into (|63p . the optimal ambiguity domain kernel (|58p is achieved. 
Validity of Eq (|59j) : Inserting the following relations into (|50p . leads to (|59p . 

Jo 



M{C Xj C* Xk )(i2-Ke)= [ 
Jo 

(Qj®Qt)(r)- 



OO c-+cu 



e c j + c fc 



U 
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